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We present a geometrical representation of sum frequency generation process in the undepleted 
pump approximation. The analogy of such dynamics with the known optical Bloch equations is 
discussed. We use this analogy to present a novel technique for the achievement of both high 
efficiency and large bandwidth in a sum frequency conversion processes using adiabatic inversion 
scheme, adapted from NMR and light-matter interaction. The adiabatic constraints are derived in 
this context. Last, this adiabatic frequency conversion scheme is realized experimentally by a proper 
design of adiabatic aperiodically poled KTP device, using quasi phased matched method. In the 
experiments we achieved high efficiency signal to idler conversion over a bandwidth of 140nm. 
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The generation of tunable frequency optical radiation 
typically relies on nonlinear frequency conversion in crys- 
tals. In this process, light of two frequencies is mixed 
in a nonlinear crystal, resulting in the generation of a 
third color with their sum or difference frequency. These 
three- wave mixing processes, also known as frequency 
up-conversion or frequency down-conversion are typically 
very sensitive to the incoming frequencies, owing to lack 
of phase matching of the propagating waves. Thus, an- 
gle, temperature or other tuning mechanisms are needed 
to support efficient frequency conversion. This difficulty 
is of particular importance when trying to efficiently con- 
vert broadband optical signals, since simultaneous phase 
matching of a broad frequency range is difficult. 

Solving the general form of the wave equations govern- 
ing three wave mixing processes in nonlinear process is 
not an easy task. Under certain conditions these can be 
simplified to three nonlinear coupled equations. Further 
simplification can be applied when one incoming wave 
(termed pump) is much stronger than the other two. In 
the "undepleted pump" approximation, two linear cou- 
pled equation are obtained rather than three nonlinear 
ones [IJ. In the case of sum frequency generation (SFG) 
process, this simplified system possesses SU(2) symme- 
try, sharing its dynamical behavior with other two states 
systems, such as nuclear magnetic resonance (NMR) or 
the interaction of coherent light with a two-level atom. 
In this letter we explore the dynamical symmetry of SFG 
process in analogy with the well known two level system 
dynamics We also apply a geometrical visualization 
using the approach presented by Bloch [3J and Feynman 
et al. [4] in NMR and light-matter interaction, respec- 
tively. The simple vector form of the coupling equation 
allows for new physical insight into the problem of fre- 
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quency conversion, enabling a more intuitive understand- 
ing of the effects of spatially varying coupling and phase 
mismatch. The utility of this approach is demonstrated 
by introducing a robust, highly efficient broadband color 
conversion scheme, based on an equivalent mechanism 
for achieving full population inversion in atom-photon 
interactions known as Rapid Adiabatic Passage (RAP) 
scheme The demonstration is experimentally realized 
using an adiabatically varying aperiodically poled KTP 
(APPKTP) in a quasi phased matched (QPM) design. 

Aperiodically poled structure have already been intro- 
duced for improving the bandwidth response of frequency 
conversion, but at a cost of a significantly reduced effi- 
ciency ini 13 HI- The broad bandwidth response is in 
particular importance for frequency conversion of ultra- 
short pulses. Chirped QPM gratings have been utilized 
to manipulate short pulses both in second harmonic gen- 
eration (SHG) [HI Uni E], difference frequency genera- 
tion (DFG) [12] , and in parametric amplification [10_, 13] . 
Recent application showed that by using a disordered 
material (random QPM), an extreme broad bandwidth 
was obtained, although at a price of severe reduction of 
the conversion efficiencies [14]. While in standard fre- 
quency conversion processes, the crystal parameters and 
the pump intensity should be precisely matched in order 
to reach high conversion efficiency, here, by utilizing adi- 
abatic frequency conversion one can still reach near to 
100% conversion efficiency over a broad wavelength and 
temperature range. In our demonstration we succeed in 
achieving near complete conversion while maintaining ex- 
tremely broad bandwidth of over 140nm. 

Let us first consider the geometrical representation of 
SFG. In the undepleted pump approximation, the pump 
amplitude is assumed constant along the nonlinear crys- 
tal, and the following normalized coupled equations for 
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FIG. 1: Bloch sphere geometrical representation of SFG in 
the undepeleted pump approximation. The amphtude of iJ2 
(bold red arrow) is much larger than the amplitudes of uji, U13, 
and therefore remains constant along the propagation. The 
energy transfer between the incoming weak field (uji) and the 
converted SFG signal(ti;3) along the propagation is dictated 
by the coupling coefficient k and the phase mismatch Ak. a) 
A periodically poled quasi phased crystal designed to achieve 
a perfect phase matched interaction along the propagation, 
(b) Geometrical visualization of the SFG dynamics on a SFG 
Bloch Sphere. Two trajectories are plotted: perfect phase 
matching (blue, torque vector points to the equator) which 
can result in efficient conversion and a constant nonzero phase 
mismatch (red, torque vector points to a point in the south 
hemisphere) , always resulting in an inefficient conversion pro- 
cess, (c) The projection of the trajectory onto the z axis yields 
the conversion efficiency. 



the signal and idler can be constructed [T]: 
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where Afc — ki + k2 — is the phase mismatch, z is the 
position along the propagation axis, k is the coupling co- 
efficient, and Ai and A3 are the normalized signal and 
idler amplitudes. The coupling coefficient and the nor- 
malized amplitudes are written below using the following 
quantities: 
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where uji and W3 are the frequencies of the signal and 
idler, respectively, ki and fcs are their associated wave 
numbers, cis the speed of light in vacuum, Ai, A2, A3 are 
the signal, pump and idler amplitudes respectively, and 
x'^) is the 2"'' order susceptibility of the crystal (assumed 
to be frequency independent). 

The dynamics of this system is dictated by the values of 
Afc and k. Full energy transfer from signal to idler (SFG 
process) or vice versa (DFG process) is usually achieved 



by demanding perfect phase matching (Afc(z) = 0), and 
also by satisfying a constraint relating the propagation 
length in the crystal with the pump intensity, such that 
K ■ z = (2n + 1) TT, where n is an integer. This is the only 
solution for a full conversion with a z-independent phase 
mismatch. 

Equations [T^ & [T]d have the same form as those de- 
scribing the dynamics in two level systems. We adopt 
Feynman et al. |1] formulation and write the dynamics 
of this problem as a real three dimensional vector equa- 
tion, which can be visualized geometrically on a sphere, 
known as Bloch sphere. In this context, we define a three 
dimensional state vector psFG — {U,V, W) as follows: 



UsFG = AIA1+AIA3, 
VsFG = i{AlAi - AIA3) 

WsFG = |A3p-|Ai|2. 
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This vector represents the relation between the sig- 
nal and idler along the crystal, and includes the co- 
herence between the modes along the propagation. In 
particular, its z component (Wsfg) gives information 
about the conversion efficiency. The south pole p = 
(0, 0, —1) corresponds to zero conversion (^3 = 0), while 
the north pole p — (0,0,1) corresponds to full con- 
version. In between, the conversion efficiency can be 
found as follow: rj = {Wsfg + 1) /2. The torque vector 
g — {Re{K}, Im{K}, Ak) represents the coupling between 
these two modes, which contains also the phase mismatch 
parameter. In this form, the loss-free evolution equations 
can be written as a single vector precession equation: 



dp, 
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g X PSFG 
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The analogy suggests the equivalence between the 
quantities in the standard two-level system and those in 
the frequency-conversion problem: The population of the 
ground and excited states are analogous to the magnitude 
of the signal and idler fields , respectively. Time evolu- 
tion is replaced by propagation in the longitudinal z-axis, 
and the detuning A is replaced by the phase-mismatch 
Afc value; The analogy is further detailed in table I. 

Presented in this context, the perfect phased matched 
solution for full conversion, as described above, has the 
same dynamical trajectory on the Bloch sphere surface 
as on-resonant light-matter interaction in atomic physics. 
This results in oscillatory dynamics between the two 
modes ("Rabi oscillations"). The case of an odd 7r-pulse 
is analogous to transferring the energy from ui (input 
signal) to the SFG/DFG output W3, i.e. the SFG state 
vector is rotated from the south pole (north pole) to the 
north pole (south pole). Any phase mismatch will lead 
to a dynamics similar to a detuned case in light-matter 
interactions, which exhibit faster oscillations and lower 
conversion efficiencies, as shown in figure 1. 

Due to the above, achievement of full energy trans- 
fer in an upcon version process is usually not robust, re- 
quiring several ingredients to be simultaneously satisfied. 
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TABLE I: Analogy between Sum Frequency Generation (SFG) process in the undepleted pump approximation and the dynamics 
of a two level atomic system, induced by coherent light. The middle column describes the optical Bloch sphere realization of 
the two level systems (notation are taken from Allen & Eberly [2]). The right column shows the analogous parameters of the 
SFG sphere realization. 
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FIG. 2: Adiabatic conversion scheme of SFG in the unde- 
pleted pump approximation, a) Continuous adiabatic varia- 
tion of the phase mismatch parameter is required. This can 
be achieved by slowly changing the poling periodicity along 
the propagation direction, b) The adiabatic following trajec- 
tory, where the torque vector initially points to the vicinity of 
the south pole, and ends up at the vicinity of the north pole, 
c) The projection of the trajectory onto the z axis yields the 
conversion efficiency. In this trajectory, phase matching con- 
dition is fulfilled at z=lcm. 



From the above analogy we can, however, adopt the RAP 
mechanism [5j, where a strong chirped excitation pulse 
scans slowly through the resonance to achieve robust full 
inversion, to the realm of frequency conversion. Thus, in 
order to transfer from Aii^z) to ^3(2) for a broadband of 
wavelengths, the phase mismatch parameter, Afc (2) has 
to change adiabatically from a big negative value to a 
big positive value (or vice versa). The adiabaticity con- 
ditions require that: 
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The first two conditions deal with the maximal magni- 
tude of the phase mismatch, which should be very large 



in comparison to the value of k, and the interaction has 
to start with a large negative (or positive) value of phase 
mismatch, and end with a large positive (or negative) 
value. The third condition is the most acute one, dealing 
with the rate of change of the phase mismatch during 
propagation. To satisfy this condition, the phase mis- 
match has to be vary slowly as compared to the square 
of the coupling term in this nonlinear process. If the 
rate of variation is not slow enough (left hand side value 
of inequality [5J; is not small enough), or the coupling 
coefficient is not large enough, this inequality won't be 
satisfied and the conversion efficiency will be poor. Fig- 
ure |2] demonstrate the case when all those constraints are 
satisfied and full frequency conversion is achieved. 

In the experimental realization of this adiabatic con- 
version scheme we choose to utilize the QPM technique 
[15'. By tuning the spatial structure of the domains, this 
technique allows not only to achieve phase mismatch that 
is close to zero along the propagation axis, but almost any 
desired function of the phase mismatched parameter. By 
proper design of the periodicity of the poling, one can 
achieve an effective phase mismatch parameter which is 
the summation of the process phase mismatch and the 
artificial phase-mismatch Ak^ (z): 



Akeff (z) 



lal 



idle 



AkA (z) . (6) 



The desired value of the phase-mismatch parameter is 
achieved by poling the QPM structure using the approx- 
imate relation: Ak\ (z) = , where A (z) is the local 
poling period. 

We demonstrated our idea by an adiabatic APPKTP 
which was designed to satisfy the constraints posed by 
equations |5] The periodicity was varied from 14.6 
to 16.2 /im along a crystal length of L = 20mm, to in- 
duce a linear adiabatic change in kf^ff(z). The design 
was carried out by a numerical simulations of the process 
propagation using the finite difference method. We as- 
sumed plain wave approximation and took x*-^-* = 32^. 
We used an OPO (Ekspla NT342) as our sources which 
produced a strong pump at 1064 nm (6 ns, 130 /iJ), and 
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FIG. 3: Conversion efficiency as a function of pump inten- 
sity with a signal wavelength of 1530nm at room temperature 
(21°C). A good correspondence between the experimental re- 
sults (shown in red) and the simulation of the design (dashed- 
dotted blue) . The maximal conversion efficiency achieved was 
74% ± 3%. 
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FIG. 4: Conversion efficiency as a function of input wave- 
length using the adiabatic APPKTP design at a pump in- 
tensity of 15 MW/crri^. A good correspondence between the 
experimental results (shown in solid red) and the simulation of 
the design (dashed-dotted blue) is shown. The low efficiency 
around 1485nm is associated with a manufacturing defect. 

a tunable signal varied from 1400 nm to 1700 nm (5 ns, 
1 nJ). The pump and the signal, both polarized in the 
extra-ordinary axis were spatially overlapped and focused 
coUinearly into the crystal having waists of ISO/im^ and 
120/im^, respectively. These values guarantee that the 
Rayleigh range is larger than the crystal length and thus 
the plane wave approximation of our simulations holds. 
We collected both the input wavelength (signal) and out- 
put SFG wavelength (idler) after their propagation in the 
crystal and recorded by an InGaAs detector, and a cooled 
CCD spectrometer, respectively. 

Our demonstration consists of two sets of experiments. 
First, to demonstrate the crucial requirement of the 
strong pump, we examined the dependence of the conver- 
sion efficiency on the pump intensity. In the second ex- 
periment we measure the conversion efficiency as a func- 
tion of signal wavelength at a constant pump intensity. 



In both cases, the conversion efficiency was measured by 
comparing the signal intensity with and without the pres- 
ence of the pump. This was checked to be completely 
correlated with the observed SFG intensity and free of 
thermal effects. 

In Fig. [3] we plot the conversion efficiency as a func- 
tion of the pump peak intensity for a fixed signal fre- 
quency 1530 nm. A very good correspondence is obtained 
with the simulation. The maximal efficiency which was 
achieved with our maximal pump intensity was 74% ±3%. 

In a second experiment, when the input wavelength 
was varied, we chose to work at a moderate pump in- 
tensity of 15 MW/cm^ and at room temperature. We 
show in Figure [4] efficient broadband conversion of over 
140nm wide (1470nm to 1610nm). This is in a good cor- 
respondence with the numerical simulation of the design 
(dashed-dotted blue), except for a small region of low ef- 
ficiency around 1485nm, which is associated with a man- 
ufacturing defect, thus causing the violation of the adi- 
abaticity condition. For comparison, a conventional ef- 
ficient converter (a periodically poled structure designed 
to achieve perfect phase matching), would lead to effi- 
cient broadband conversion only for 2nm bandwidth. 

Other experiments showing the robustness of the adia- 
batic conversion scheme to errors in incidence angle, crys- 
tal temperature and pump frequency (as compared to the 
periodically poled design), and exhibiting color tunabil- 
ity of the conversion band by temperature control will be 
presented in a forthcoming publication. 

In conclusion, we show that the use of a geometri- 
cal Bloch sphere visualization to describe evolution of 
the complex mode amplitudes, gives information on their 
relative phases as well as intensities. We demonstrated 
high efficiency wavelength conversion for a wide range 
of frequencies, by using the adiabatic conversion method 
adapted from other well known two level systems. It is 
important to note that quasi phase matched crystal is 
only one possible realization; the same mechanism can 
be applied, for example, by inducing a temperature gra- 
dient on a perfect phased matched crystal, or by any 
other mechanism fulfilling the adiabatic constraint of Eq. 
|5] The present scheme can be utilized for efficient fre- 
quency conversion of broadband signals as well as ul- 
trashort pulses. This analysis holds promise for ultra- 
short pulse conversion in a wide range of frequencies from 
UV to far IR. It may be particularly useful in the effi- 
cient upconversion of weak infrared signals to the near- 
infrared or visible range, used in atmospherical observa- 
tion [TSJ [T7] . The generalization of this analogy to higher 
order nonlinear processes interaction can be made, and 
known schemes from atom-photon interaction or NMR, 
such as Stimulated Raman Adiabatic Passage (STIRAP) 
[TH] can be adopted. 
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